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Big Bang Nucleosynthesis (BBN) has increasingly become the tool of choice for investigating the 
permitted variation of fundamental constants during the earliest epochs of the Universe. Here we 
present a BBN calculation that has been modified to permit changes in the QCD scale, Aqcd- The 
primary effects of changing the QCD scale upon BBN are through the deuteron binding energy, 
Bu, and the neutron-proton mass difference, Sm np , which both play crucial roles in determining the 
primordial abundances. In this paper, we show how a simplified BBN calculation allows us to restrict 
the nuclear data we need to just Bu and 8rn np yet still gives useful results so that any variation 
in Aqcd may be constrained via the corresponding shifts in Bu and Sm np by using the current 
estimates of the primordial deuterium abundance and helium mass fraction. The simplification 
predicts the helium-4 and deuterium abundances to within 1% and 50% respectively when compared 
with the results of a standard BBN code. But Aqcd also affects much of remaining required nuclear 
input so this method introduces a systematic error into the calculation and we find a degeneracy 
between Bu and 5m nv . We show how increased understanding of the relationship of the pion mass 
and/or Bo to other nuclear parameters, such as the binding energy of tritium and the cross section 
of T + D He + n, would yield constraints upon any change in Bu and 5m np at the 10% level. 

PACS numbers: 99.99 



I. INTRODUCTION 



Speculation that fundamental constants may vary with 
time began as early as the 1930's pj and, although there 
is no mechanism for such time variation in the context 
of the standard model of particle physics, recent obser- 
vations have created renewed interest in this idea. Webb 
et al. report observations of quasar absorption lines 
at redshift of z=l-2 that suggest the fine structure con- 
stant, a, may have been smaller at this time. The effect 
is at the level of one part in 10 5 and further analysis of 
a new data set by the same group gives similar results 
although Bahcall, Steinhardt and Schlegel Q using a 
different analysis method find a different limit. 

Independently of these observations however, it is in- 
teresting to consider what happens if the fundamental 
constants were different at earlier time than they are to- 
day. There are many suggestions for beyond the standard 
model theories that could accommodate time variation in 
the fundamental constants and link the changes in some 
constants to others, e.g. Langacker [3, though in general 
studies have typically derived constraints on one constant 
when all the others are fixed. If these constants were dif- 
ferent at an early epoch than they are today, their relative 
shifts can only be determined by the underlying theory 
which causes the changes. In the context of grand unifies 
theories, all couplings are correlated via the GUT scale 
parameter, which varies with time, and relative ratios 
are determined by the renormalization group equations. 
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The relationship between the shift in the fine structure 
constant and the vacuum expectation value of the Higgs 
boson is very model dependent even in the context of 
grand unified theories. For some supersymmetric theories 
it is as large as Av/v ~ 70Aa/a. While the limits on 
the variation of the fine structure constant derived from 
quasar spectra (which were formed ~ 1 billion years after 
the Big Bang) are severe and would seem to limit the 
amount we could expect in a variation of Aqcd at that 
epoch (for example from equation ^) they by no means 
rule out the possibility of much larger changes in Aqcd, 
a or indeed any other constant at the much earlier phase 
of BBN unless a model is provided that allows one to 
calculate how their value at one epoch is related to the 
value at another. 

To constrain the permitted variations one can look 
to several places to confront theory with data and the 
strongest constraints upon the time variation of funda- 
mental constants are expected to emerge when the ob- 
servables come from events in the distant past. Such 
places include: the Oklo nuclear reactor, quasar absorp- 
tion spectra, the Cosmic Microwave Background and Big 
Bang nucleosynthesis. There have been a number of stud- 
ies that consider the variation of particular fundamental 
constants in these scenarios. Data from The Cosmic Mi- 
crowave Background anisotropy measurements and BBN 
were used by Yoo and Scherrer Q for the Higgs vacuum 
expectation value while the effect of changing the fine 
structure on Big Bang Nucleosynthesis has been exten- 
sively studied by Bergstrom, Iguri and Rubinstein H and 
Nollett and Lopez Flambaum and Shuryak fla, IT^j 
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give constraints upon the quark masses and Aqcd af- 
ter using BBN and limits derived from the Oklo nuclear 
reactor have also been studied 0, 0] . In a study that 
permitted changes in all the gauge and Yukawa couplings 
by relating them to the evolution of a single scalar field 
(the 'dilaton') Ichikawa and Kawasaki again use BBN 
to limit the variations of these couplings. Fundamental 
coupling constraints on BBN were also studied by |l2j. 

In this paper we revisit the study of the influence of 
the strong coupling on Big Bang Nucleosynthesis yields. 
BBN occurs in the first few minutes and hours after the 
Big Bang and hence we would naively expect the maxi- 
mal difference from the current values of the fundamental 
constants if indeed they vary with time. A variation in 
h-QCD is much more difficult to implement than changes 
in most other constants because of the uncertainty of how 
the difference would impact the nuclear data required to 
make a prediction. Nevertheless, that is our intention in 
this paper. In the next section we show that of all the 
inputs to the calculation, it is the deuteron binding en- 
ergy, _Bd> and the neutron-proton mass difference, 5m np 
that play the most crucial role in determining the primor- 
dial helium and deuterium abundances - the two nuclei 
with, currently, the two most accurately known primor- 
dial abundances with which to compare the calculation. 
In the following section we discuss how the variation in 
fundamental parameters is related to the deuteron bind- 
ing energy and the neutron-proton mass difference and 
discuss how these quantities enter into the important re- 
action rates. In the fourth section we describe a calcula- 
tion of Big Bang nucleosynthesis with a modified reaction 
network and show explicitly the effects of 5m np and -Bd 
on abundance yields before deriving constraints on these 
quantities. In the last section we give our conclusions 
and point to where further understanding of the interac- 
tions/structure of nuclei can improve our results. 



II. STANDARD BBN 



characteristics of the setting: the duration, the energy 
available and the interactions between the nuclei. Sim- 
plistically the primordial abundances are determined in 
three (somewhat) distinct phases whose boundaries are 
determined by just two nuclear parameters: the neutron- 
proton mass difference 5m np and the binding energy of 
the deuteron Bd. In each, the behavior of the nuclei is 
distinguished by their interactions and, in order to set the 
stage for later discussions, we shall skip briefly through 
each pointing out the importance of these two parame- 
ters. 



A. 8rn np : Weak Equilibrium 

Throughout the entire evolution of the baryons in BBN 
it is the relativistic particles that drive the expansion, 
provide the thermal bath in which the nuclear reactions 
take place, and set the time-temperature relationship. 
At a sufficiently early epoch all particles are in ther- 
mal contact through electromagnetic and/or weak inter- 
actions and so possess a common temperature T. At a 
temperature of 10 MeV, a common initial temperature 
for BBN calculations, the radiation includes: photons, 
electron/positrons and three light neutrinos. The in- 
consequence of the nucleons / nuclei for cosmology at this 
epoch is implied by the smallness of the ratio of baryon 
and photon number densities, denoted by 77, which is 
V = ns/n-y ~ O(10~ 10 ). Any chemical potential of 
the photons is driven to zero from such rapid processes 
as double Compton scattering, in the absence of any 
(significant, non-thermal) source of neutrinos the ratio 
S,i = fJ-i/Ti is constant, and the chemical potential of the 
electrons/positrons is set by the proton density and is 
therefore extremely small, £ e ~ 77 [Ty| . The energy den- 
sity of the Universe is thus determined from the common 
temperature and the chemical potential of each relativis- 
tic fluid component so that the expansion rate of the 
Universe, denoted by the Hubble parameter H , is simply 



BBN represents the marriage of nuclear physics with 
cosmology and, in comparison to the majority of nucle- 
osynthesis settings, the calculation is relatively simple 
since there are no production zones to deal with (the 
whole Universe participates) and no spatial gradients of 
any kind that lead to a transport of entropy, momentum 
or mass. But the continual dilution and cooling of the 
cosmic fluid certainly does not mean the Universe is in a 
steady-state during these early phases. The expansion is 
driven by the energy density of the relativistic particles, 
the nucleons/nuclei represent only a small fraction of the 
total, and the inexorable decrease in temperature and 
density means that the nuclear reactions can occur for 
only a brief period. Despite their cosmological inconse- 
quence the (inferred) abundances of the nuclei at the end 
of BBN represent the best set of observables that probe 
sub-horizon scale physics at this early epoch. Like all nu- 
cleosynthesis mechanisms BBN is sensitive to three key 
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and the age of the Universe is related to the temperature 
through 



tT 2 = 0.74 MeV 2 s, . 



(3) 



when all neutrino chemical potentials are zero. 

The neutrons and protons are also held in chemical 
equilibrium via such weak interactions as 



n + e p + v e 



(4) 



From the requirement fi n — /_t e = /_t p — fi Ve the neutron to 
proton ratio F is 
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As the Universe cools the plethora of equilibria estab- 
lished above 10 MeV is broken by the increasingly in- 
frequent weak interactions. The simplicity of a single 
temperature for every fluid constituent falters when the 
neutrinos 'decouple'. The ebb of the weak interactions 
occurs when the interaction timescale becomes longer 
than the age of the Universe so that the neutrino de- 
coupling temperature represents the point at which the 
two are equal. Below this temperature a neutrino is, 
on average, unlikely to ever experience another interac- 
tion that would allow energy to be transferred from one 
fluid component to another. For v p and v T the decou- 
pling temperature is T ~ 3.3 MeV, for v e it is slightly 
lower at T ~ 2 MeV because v e may also interact with 
the electron/positron fluid via exchange whereas the 
other two flavors do not. Despite the fact that the three 
neutrinos no longer interact with any other fluid compo- 
nent there is no change in the evolution of the Universe. 
Since there was no change in the number of degrees of 
freedom in any of the fluid components during decou- 
pling the three flavors all possess the same temperature 
T v equal to the electromagnetic temperature T~. Neu- 
trino decoupling has no direct impact upon the nucleons 
at this point in their evolution, it's importance is resur- 
rected during a later epoch. 

Concurrent with neutrino decoupling the nucleons also 
undergo an equilibrium crisis. The neutron/proton in- 
terconversion reaction n <-> p, which actually represents 
three processes, 



n «-> p + e + v e 
n + e «-> p + v e 
n + v e «-> p + e, 



(6a) 
(6b) 
(6c) 



is also governed by the weak interaction and is the only 
process that can significantly alter the neutron/proton 
ratio. 

At temperatures of ~ 10 MeV the three body reaction 
p + e + D e — ► n and its inverse, neutron decay, both oc- 
cur at a rate smaller than the Hubble parameter and are 
therefore incapable of affecting the neutron to proton ra- 
tio. In contrast the two-body reactions (j6b|) and (|6*c|l are 
sufficiently rapid initially to establish an equilibrium but 
below T 7 ~ 1.2 MeV, a temperature we shall denote by 
T n p, the rates become smaller than the Hubble rate and 
are therefore incapable of maintaining this equilibrium. 
The neutron to proton ratio is said to 'freeze out' though, 
in truth, even if all complex nuclei were prohibited and 
neutrons did not decay the ratio would not become a 
constant until much later [20, 21]. 

The departure from neutron/proton weak equilibrium 
sets the first milestone in the path to BBN. The bound- 
ary between the equilibrium phases prior to ~ 1 MeV 
and the non-equilibrium phase that follows is determined 
essentially by the single parameter Sm np and this quan- 
tity also determines the freeze-out ratio of neutrons and 
protons. The freeze-out ratio sets the upper limit to the 
number of neutrons that can participate in BBN and will 



therefore limit the primordial abundance of every com- 
plex nucleus. As the Universe continues to evolve below 
~ 1 MeV the n/p ratio decreases from the freeze-out ratio 
because of both the fading residual two-body interactions 
and the emerging importance of neutron decay: the ex- 
tent to which this limits the number of neutrons that will 
go on to participate in BBN proper is determined by the 
second parameter By>. 



B. B u : NSE 

There are several key events that occur during this 
second stage of the maturing Universe. The first is to 
the background fluid when the number/energy density 
of e, e begins to depart from the relativistic value. This 
occurs at roughly T 1 ~ 800 keV and does not cease until 
T 7 ~ 10 keV when the annihilation rate becomes smaller 
than the Hubble rate. The electron/positron annihila- 
tion deposits energy and entropy only into the photon 
gas, the neutrinos have decoupled and cannot share in 
this energy release. The change in the number of electro- 
magnetic degrees of freedom leads to an increase of T 1 
relative to T„ in order to maintain the entropy within the 
co-moving volume though the increase is never capable 
of reversing the redshifting due to the expansion. After 
annihilation the electromagnetic temperature T 7 is re- 
lated to the neutrino temperature T„ via the well known 
T 3 /T 3 = 4/11 and the time-electromagnetic temperature 
relationship evolves to become 



tT 2 = 1.32 MeV 2 s 



(7) 



again with the standard assumption of zero neutrino de- 
generacy. The changing time-temperature relationship is 
important here because it determines how much neutron 
decay can occur before nuclei begin to form. Hereafter we 
drop the subscript for T 7 so that whenever we mention 
temperature it is always the electromagnetic. 

So far we have made no mention of the complex nuclei. 
The temperatures are so high that their abundances are 
suppressed relative to the free nucleons but of course the 
nuclear reactions that form them do occur. In contrast 
with the n <-> p interconversion processes the nuclear re- 
actions, such as n + p <-> D + 7, are rapid at T np and so 
the abundances reach, and maintain, chemical/nuclear 
statistical equilibrium (NSE). In equilibrium the abun- 
dance 1 , Ya — ua/tlb, of the complex nuclei A is derived 
from fiA = Z jU p + (A — Z) /x n so inserting the expressions 
for the Boltzmann number density we find. 
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1 The term 'abundance' is also used for the ratio Ya/Yh 
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From this equation we see that for a temperature of T ~ 
1 MeV the abundance of deuterons is Yd ~ 10 -12 . We 
can rewrite equation (JHJ to illustrate what is happening in 
this temperature regime, by replacing the thermal factors 
with the Yd so that 
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This equation makes it much clearer that the abundance 
of nucleus with mass A + 1 relative to nucleus with mass 
A is smaller by approximately Yd indicating just how 
severe the suppression is when Y D is small. 

There are many different reactions in which the nuclei 
participate. For BBN the most important are [22| 
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D + D <-> n+ 3 He 
D + D <-> p + T 
D + p <-> 7 + 3 He 



T + D 

T+ 4 He 



n+ 4 He 
7+ 7 Li 



J He 
3 He- 
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4 He 



3 He+ 4 He <-> 7 + 7 Be 



7 Li + p <-> 2 4 He 
7 Be + n <-> p 



7 Li 



(10) 



(11a) 
(lib) 
(11c) 



(12a) 
(12b) 



(13a) 
(13b) 
(13c) 

(14) 
(15) 



In every case the rate of change of abundance of a 
nucleus j, a member 2 of the set of products {P}, from 
any given reaction involving the reactant set {R} is 



1 dYj Yp 



(16) 
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where the u's are the stoichiometric coefficients and 
Tr—>p is the rate of the reaction per unit abundance of 
the reactants. In equilibrium, the production and de- 
struction of the nuclei by any one reaction are almost 
equal with only a small difference between them. This 



allows us to relate the forward, rp_>p, and back reaction, 
Tp^ji, rate coefficients i.e. 
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If we insert into this expression the Bolztmann distribu- 
tion of the number densities then we find 

Tp^ R oc Tp^p exp [ 118) 



T 



r / Q 

oc Tp^p exp ( — 



(19) 



where Mr and Mp (/zp and /xpjare the sum of the reac- 
tant or product particle masses (chemical potentials) and 
Q is the difference. If the reaction involves only nucle- 
ons/nuceli and photons then /ip = fip and the l Q value' 
is simply the difference in the total binding energy of 
reactants and products, but in those cases where the re- 
action involves particles we have not included in {R} and 
{P} the equality is not ensured. An analytic discussion 
of the effects BBN from the point of view of the reaction 
rates can be found in [23| . 

As the Universe cools the reactions are unable to pro- 
cess the nuclei leading to a departure of the abundance 
of each from the NSE abundance: the heavier nuclei de- 
parting earlier than the lighter. The ratio of the NSE 
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FIG. 1: The ratio Y i NSE /Y i , the NSE abundance to the actual 
abundance, of D, T, 3 He and 4 He as a function of the photon 
temperature when the baryon/photon ratio is 5.6 x 1(F 10 . 



abundance and the actual abundance from a standard 
BBN code is shown in Fig. QJ and was first discussed in 
Smith, Kawano & Malaney [22j. The fi gure shows that 
the 4 He abundance departs from NSE when T ~ 600 keV 
while 3 He and T depart at T ~ 200 keV. Below 
T ~ 200 keV only the abundance of the deuteron is given 
by its NSE value, the rest are many orders of magnitude 
smaller. The departure from NSE is not because the rates 
fall beneath the Hubble rate but rather the departure oc- 
curs because the amount that can be provided from the 
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reactions falls short of the amount required to remain in 
equilibrium. For any particular nuclear species the for- 
mer is simply the sum of all the relevant production and 
destruction reactions listed in (|1(J|) through (|15|l while the 
latter can be found from equation J3J) (since deuterons are 
only beginning to form, we ignore the change in Y p and 
Y n and remembering ub oc T 3 ), 



dY A 



_ Y_A_ dT_ 

dt T dt 



3(^4-1) 



Ba 
T 



(20) 



The departure from NSE has a major impact upon the 
reactions. The rate of change of nucleus k in, for example, 
a two-body reaction i + j ^ k + I with i =/= j, k I, is 



dYk 
dt 
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(21) 
(22) 



with the superscript NSE indicating the equilibrium 
abundance, is now very lopsided because Yk/Y l ^ fSE <C 1. 
Essentially there is no destruction of k via this reaction, 
the back reaction has switched off. 

By T ~ 200 keV the only compound nucleus in NSE is 
the deuteron. Its abundance is rising rapidly and BBN is 
said to begin when the nuclear reactions finally begin to 
process the deuterium leading to the final departure from 
NSE. The increasing significance of such reactions as D + 
D — > T + p halts the rise in the deuterium abundance 
and Bernstein, Brown & Feinberg define the temperature 
at which BBN begins, Tbbn, as the point where the 
deuteron abundance reaches its peak i.e. dY^/dt = 0. 
Using this definition they find Tbbn ~ -Bd/26 = 86 keV, 
or t'BBN = 180 s. 

The second stage to BBN, from 1.2 MeV to 86 keV, 
is characterized by departure from equilibria. At its in- 
ception non-equilibrium is limited to the neutron/proton 
ratio: by its end all other nuclei have fallen out of equilib- 
rium with deuterium the last to succumb. Formation of 
the heavier nuclei in significant amounts begins to occur 
towards the end of this stage and indeed the abundance of 
4 He can already be larger than D by Tbbn ■ What Tbbn 
represents is the point where the (leaky) dam bursts so to 
speak. The initiation of the next stage of BBN proper is 
controlled by the deuteron binding energy £?d and there- 
fore, in conjunction with 5m np , the initial conditions for 
BBN proper are essentially controlled by just these two 
parameters. 



C. BBN proper 

From the departure of n/p weak equilibrium at T ~ 
1 MeV to the inception of BBN at T ~ 100 keV some 
of the free neutrons have decayed and reduced the pool 
available to be assimilated. But once the balance has 
been tipped in favor of complex nuclei the free nucleons 
are rapidly dragged into and through A = 2 and neutrons 



become stabilized. This phase is BBN proper and dur- 
ing it the abundances of the complex nuclei can become 
very large as the nuclear reactions process them. But by 
the time the Universe has cooled to the point where the 
reactions have ceased virtually all the neutrons present 
at Tbbn now reside in helium-4 with a small fraction 
in the trace abundances of the other nuclei. The trace 
abundances of the intermediary nuclei D, T and 3 He are 
the ashes of processing from free nucleons to 4 He . Be- 
cause helium-4 essentially acts as the end-point of the 
reactions its abundance is very insensitive to the exact 
details of the nuclear reactions that lead to its forma- 
tion. The abundance of 4 He instead probes the much 
earlier epoch of n/p freeze-out and the time delay until 
BBN commences. In contrast, the final abundance of the 
intermediary nuclei D, Tand 3 He are strong functions 
of the nuclear physics. The efficiency of the processing 
from free nucleons to 4 He depends crucially on the con- 
ditions during BBN, the temperature and the interaction 
cross sections, so their abundances can vary markedly if 
BBN begins at a higher temperature or the cross sections 
change. 



III. NUCLEAR REACTION RATES AND A QCD 

As discussed above, the deuteron binding energy and 
the neutron-proton mass difference are two of the most 
important parameters that control Big Bang Nucleosyn- 
thesis. Each of these depends on the QCD scale, although 
each in a different way. 

Shifts in the QCD scale might not, at first glance, be 
expected to have significant effects. This is because most 
nuclear quantities, e.g. nuclear masses, nucleon masses, 
will shift together with this scale and therefore not be 
observable. However, there is no reason to expect that 
all the fundamental constants shift at the same rate, and 
the relative change may cause observable effects. For 
example, while the masses of the proton and neutron 
are mainly determined by the QCD scale, the difference 
between the neutron and proton masses is dominated by 
electroweak effects. As a consequence, the ratio (m n — 
m p )/m n may be sensitive to changes in the fundamental 
couplings. Furthermore the pion mass will not shift in the 
same way as the nucleon or nuclear masses, and the pion 
mass is an important ingredient in setting the deuteron 
binding energy. 

The neutron-proton mass difference is approximately 
given by pi} 



8m 



np 



M n - M p = m d - m, 
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The coefficient of a in the electromagnetic contribution, 
M e i m , is determined by strong interactions and therefore 
proportional to Aqcd- We use the estimate aM e i m ~ 
0.76 MeV, while the difference in the down and up quark 
masses is approximately 2 MeV [24J. Other calculations 
of aM e i m yield different results, for example [25|. Shifts 
in the mass of the up and down quark masses depend 
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on the vacuum expectation value of the Higgs and their 
respective Yukawa couplings. Specifically, 



A{8m np ) = - I 



a 



AA 
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Here the y's are the Yukawa couplings and v is the vac- 
uum expectation value (vev) of the Higgs boson. From 
this equation, it can be seen that the shift in the neutron- 
proton mass difference is determined by several funda- 
mental parameters, not just Aqcd- 

In principle, we should define a scale which remains 
fixed, while others vary. This could be a scale such as 
the GUT scale. In our calculations in the next section 
however, we only vary Aqcd so this scale could even 
be the electron mass. A discussion of scales as well as a 
preliminary estimate of the deuteron binding energy may 
be found in [2(j 

The scale Aqcd also has great importance in deter- 
mining the binding energy of the deuteron. The small 
binding energy prevents a significant abundance of deu- 
terium from building up early in big bang nucleosyn- 
thesis. And therefore, as discussed above, the deuteron 
binding energy controls the nuclear flow to elements with 
greater mass number. The long-range part of the nuclear 
force is governed by one-pion exchange and therefore the 
deuteron binding energy is sensitive to pion properties. 
The mass of the pion is determined by the Gell-Mann- 
Oakes-Renner relation, 1271 



Here / w is the coupling of the pion to the axial current, 
and < qq > is the quark condensate. 

Shifts in the pion mass are then calculated in a simi- 
lar manner to Eq (|24(l and determined by the shifts the 
Yukawa couplings, v and Aqcd- Since f n oc Aqcd and 
< qq >oc A 3 QCD , 



Am, 



AAqcq Av Ay u + Ay d 
Aqcd v Vu+Vd 



(26) 



llm\ = (m u + m d ) < qq > 



(25) 



Recent studies of the dependence of the nuclear potential 
on the mass of the pion have shown a strong variation 
with the mass of the pion [28l 12^. l30| . According to 
these authors, changes in the pion mass of a few percent 
could lead to changes in the deuteron binding energy of a 
factor of two. Although a relationship between the pion 
mass and the deuteron binding energy exists, [2^, I29I 
l3fl | , several parameters in this relationship are uncertain. 
These include parameters that describe the pion nucleon 
coupling and the coefficient of a quark mass dependent 
four nucleon operator. 

Next we turn to the reaction rates themselves that are 
most important for determining the deuterium and he- 
lium abundances. 

The first is the neutron-proton interconversion reaction 
listed in (|6a|l through ifficjl . The rates for the reactions 
must include the temperature chemical potentials of the 
electrons and neutrinos which also participate in the re- 
actions. In the Born approximation the neutron/proton 
interconversion rates are 
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(27) 

/ dE e E e p e [(E e + 5m np ) 2 f s (E e ) f D {E e + 5m np ) + (E e - 5m np ) 2 f e {E e ) j v {E e - 5m np )] 

J m P 

(28) 



where f(E) is the Fermi-Dirac distribution, f(E) its com- 
pliment and the subscripts indicate the chemical poten- 
tial and temperature to be used in /. Though these two 
expressions capture the essence of how Sm np enters into 
these rates there are a number of corrections that must 
be taken into account before the n <-> p reaction rates 
reach sufficient accuracy 0, 0, 0] given the accuracy 
of the data we will eventually use in our constraints. 

The next is the reaction which creates deuterium, 
n + p — > D + 7, which is well studied in a low energy ef- 
fective field theory without pions 0, . These authors 



give an expression for the cross section as a function of 
the deuteron binding energy, the scattering length in the 
singlet channel, the phase shift, and so on. In principle 
each of these parameters should be treated as free, how- 
ever, the one which has by far the largest leverage on the 
rate is the binding energy of the deuteron. In our calcu- 
lation, we use the expression for the cross section given 
by 32] and integrate over the thermal distribution of the 
particles in the manner prescribed in Fowler, Caughlan 
and Zimmerman |33| ]. 

The other crucial reaction is the one that destroys most 
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of the deuterium, D + D — > T + p during BBN proper. 
Ideally, we would like to use the results of a nuclear ef- 
fective field theory calculation in order to determine the 
dependence of this rate on the pion mass, but there are 
no effective field theory calculations available yet. 



IV. NON-STANDARD BBN 

Standard BBN is well studied problem and even 
though the compatibility of its predictions with the ob- 
served abundances remains a contentious issue much of 
the recent efforts have been towards using it as a test of 
the state of the Universe at the earliest epochs. Modi- 
fying standard BBN to include new effects can be rela- 
tively painless: for example, the use of BBN as a probe 
of the number of light neutrino species or the chemical 
potentials of the y, and/or r neutrino flavors is simply 
a case of modifying the energy density and hence the 
expansion rate with no direct influence upon the nu- 
clear physics 0, EE EH • Similarly the expansion rate 
can be modified by changing the cosmological equations 
[37l l38| . introducing extra energy density in the form 
of 'Quintessence' Hi IH E3, Ei 

and even quite general 
evolutions have been explored [4j|. The effects of new 
physics from scenarios such as a non-zero electron neu- 
trino chemical potential |H El EjMElJIlJ^Js^r 
neutrino oscillations/mass/decay [5lL 152. 153L 154 l55l l56j . 
show up primarily in the neutron-proton interconversion 
reactions. There the known expressions for these reaction 
rates permit a high degree of confidence that all the con- 
sequences of the new physics have been accurately taken 
into account. 

By far the most difficult proposals to implement 
touch every nuclear reaction. For example, the impact 
of a varying fine structure constant was examined by 
Bergstrom, Iguri and Rubinstein [5Sj. later improved 
upon by Nollett and Lopez 59], and to examine its ef- 
fects these authors had to recalculate every relevant rate. 
In this type of study one is forced to make a number 
of approximations because we do not have a complete 
understanding of how a change in a fundamental con- 
stant alters low energy nuclear physics parameters such 
as the binding energies and cross sections. An alter- 
native approach would be to map the parameter space 
from a single, or handful of, fundamental constant(s) to 
a greater number of nuclear physics parameters, and then 
constrain the nuclear physics parameters, while treating 
them as independent. This way we avoid the uncertainty 
in how changes in the fundamental constants manifest 
themselves in nuclear quantities but this comes at the 
expense of exploring larger regions of parameter space. 
Despite this trade-off, we shall turn in this direction. 

The number of nuclear physics parameters is huge but, 
as we have been at pains to stress, the most important 
two are 5m np and £?d and indeed we can limit our pa- 
rameter space to just these two quantities. Before we set 
out our plan for constraining 5m np and Bd we examine 



why we can restrict the number of nuclear parameters to 
just this limited set yet still obtain good predictions for 
the results of BBN and how these two parameters will 
affect the predicted primordial abundances. 



A. Simplifying BBN 

A prediction for the mass fraction of helium-4 is reli- 
ably calculated by simply counting the number of neu- 
trons that become stabilized inside nuclei. This approx- 
imation for the helium-4 mass fraction emerges from the 
empirical results of standard BBN discussed above where 
we showed that all the neutrons that survive BBN re- 
side in helium-4. To determine the number of stabilized 
neutrons we must first be able to calculate the neutron- 
proton ratio during the earliest, weak-equilibrium, phase 
of BBN and then, secondly, follow it through to the in- 
ception of BBN proper. The first is relatively easy to 
implement since we have known expressions to use for 
the neutron-proton interconversion rates. As we men- 
tioned earlier, there are many corrections to equations 
(|27[1 and (|28[1 that must be included in order to improve 
their accuracy. Of these we only explicitly included the 
radiative electromagnetic corrections and do not intro- 
duce the corrections for the finite mass or thermal ra- 
diative corrections. Since they are relatively small when 
compared with the effects we want to consider, we take 
them into account only by an overall scale factor normal- 
ized to the measured neutron lifetime. The second step 
means that we must be able to follow the neutrons into 
the first reaction that must occur in processing the nu- 
clei: n + p «-> D + 7. As we mentioned in section Hlfll an 
analytic expression for the cross section of this reaction is 
available from Rupak [33 ] which shows how the deuteron 
binding energy enters into this important quantity. 

As we learned from standard BBN the final abundance 
of 4 He is very insensitive to the exact rate of its formation 
so if we truncated our reaction network at this step we 
could still get a good estimate for the helium-4 mass frac- 
tion and thus we only need two parameters for the calcu- 
lation: the neutron-proton mass difference 5m np and the 
deuteron binding energy Bu- But with two parameters 
(three if we include the baryon to photon ratio 77) and 
one prediction a degeneracy is expected to arise between 
5m np and Sd in the abundance of 4 He . Therefore, it 
is worth endeavoring to make a prediction for the abun- 
dance of at least one other nucleus. The most obvious 
candidate is deuterium because its primordial abundance 
is also a function of these same two parameters. To make 
a prediction for the amount of unprocessed deuterium we 
have to follow not only how this nucleus is formed cour- 
tesy of the n + p <-> D + 7 reaction but also how much is 
destroyed. The majority of the destruction of deuterium 
occurs via the reactions (|Ila|l . l|llb(l and l|I2a|) with (|f I c|) 
and i|13b|l following behind in importance |22|. As we saw 
in equation (|I9|I the inverse rate is related to the forward 
rate by a term that is proportional to exp(— Q/T) and Q 
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may be expressed in terms of the binding energies. By 
including these reactions we are adding the tritium and 
helium-3 binding energies explicitly to the 5m np and £?d 
parameter list through this exponential dependence. (In 
addition, there may be additional dependence of these 
binding energies in the forward cross sections, although 
there is no explicit calculation currently available.) But, 
here another result from standard BBN comes to our res- 
cue: the nuclear abundances of all the nuclei heavier than 
D during BBN fall far below their nuclear statistical 
equilibrium values and the inverse reactions essentially 
switch off long before BBN begins. This allows us to use 
tritium and hclium-3 as the substitutes for helium-4 if we 
can assume that any neutron that ever forms one of these 
nuclei is prevented from ever participating in another re- 
action. If we are able to ignore all inverse rates then we 
are eliminating the explicit dependence upon the bind- 
ing energies of tritium and helium-3 (although the cross 
sections will still implicitly be functions of these quanti- 
ties). Furthermore, if we truncate our reaction network 
at these two A = 3 nuclei then we need only be concerned 
with reactions Ijllajl and (|llbjl . With the absence of an 
analytic expression to use, we take the cross sections of 
these two reactions to be to be related to the size of the 
deuteron radius i.e. a ~ \/{m^ B-q). In terminating the 
network so early we miss the significant destruction of 
deuterium from T + D — > n + 4 He and the less important 
3 He + D — > p + 4 He but nevertheless, we can use the 
abundance of D from the calculation as an estimate of 
the primordial abundance even though this approxima- 
tion is much cruder than the one used for the primordial 
4 He abundance. 

Our approximation to BBN obviously introduces un- 
certainty into its predictions but the agreement with 
standard BBN calculations is remarkably impressive. As 
a function of the baryon to photon ratio rj the helium-4 
abundance is systematically too low by < 1% compared 
to a standard calculation. The reason lies in the fact 
mentioned previously: the abundance of 4 He surpasses 
the abundance of D at Tbbn, the beginning of BBN 
proper, but since we have essentially removed this nu- 
cleus from the network we overestimate the amount of 
neutron decay and therefore calculate a lower than ex- 
pected yield of helium-4. In contrast deuterium is sys- 
tematically overestimated but this time by ~ 50%. Ne- 
glecting T + D — > n + 4 He has introduced insufficient 
deuterium destruction , which is of order the observation 
uncertainty. We will show later how powerful the use of 
deuterium could be if we had more information about 
how other nuclear cross sections/binding energies vary 
with the pion mass. 



B. The effect of changing Bd 

The effects of a change in B-q occur through both the 
change in the NSE deuteron abundance and the cross 
sections that process D to heavier nuclei. As B-q in- 



creases the deuteron becomes more stable and conse- 
quently more difficult to disassociate. Therefore we ex- 
pect as Bd increases BBN will commence at a higher 
temperature/earlier epoch. Perhaps counterintuitively 
the increased stability of the deuteron leads to a decrease 
in the primordial deuterium abundance. The rate coef- 
ficients are functions of both the temperature and the 
powers of the Baryon density tib so if Tbbn increases 
interactions occur more rapidly even if the cross sections 
have changed. This leads to a more efficient processing of 
the neutrons to 4 He so that the amount of unprocessed, 
and thus the primordial, deuterium decreases. From a 
more efficient processing of the nucleons we would expect 
an increase in the primordial 4 He abundance and this is 
enhanced by the shorter interval between n/p freeze-out 
and the beginning of BBN and consequently less n-decay. 
These expectations are confirmed in figure J2J) where we 
plot the deuterium abundance Yd and helium-4 mass 
fraction, confusingly denoted also by Y, as functions of 
the temperature for ABd /Bd = and ABd/Bd = 0.3. 
At temperatures above 100 keV when deuterium is in 

10-2 I ' ' ' ' ' ' 1 
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T(eV) 

FIG. 2: The deuterium abundance, Yd, top panel, and the 
helium-4 mass fraction, Y, bottom panel, as functions of the 
photon temperature for ABd /Bd = and ABd/Bd = 0.3 
when the baryon/photon ratio is 6.14 x 10 -10 . 



NSE the offset due to the change in its binding energy 
is clear, as is the shift in the peak deuterium abundance, 
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and hence Tbb n ■ The lower figure shows that the change 
in the deuterium binding energy leads to an increase in 
the final helium-4 mass fraction as predicted. 

C. The effect of changing Sm np 

There are two effects that become apparent when we 
change the neutron-proton mass difference: an altered 
neutron lifetime and a change in the n/p freeze-out ra- 
tio. Before we add the radiative corrections, the neutron 
lifetime r n is simply 

1_ ^ G\ (C£ + 3Cp 
r n ~ 2tt 3 

x / dE e E ePe (E e -Sm np ) 2 (29) 

= Sm^ p I(m e /Sm np ). (30) 

Since function I(m e /5m np ) varies less rapidly than <5m^ p 
as 6m np increases the lifetime drops considerably. We 

3 Although Bernstein, Brown & Feinberg 20] ignore n-decay in 
deriving their expressions they still require neutrons to decay 

I 



where the function f(m e /Sm np ) varies slowly at 
m e /Sm np w 1/2. The change of Y* in i|31|) is there- 
fore primarily due to the effect of the r Sm^p factor in 
the denominator of the exponential. From equation l|29() 
we know that the lifetime scales as ~ 1/ <5m^ p so the fac- 
tor in the exponent is proportional to l/<5m„„. As 5m n p 
increases the exponential term in i|31|) fades more rapidly 
yielding a faster convergence with q and so confirming our 
expectation of a lower freeze-out abundance of neutrons. 

Finally, as a consequence of the lower freeze-out neu- 
tron abundance the NSE deuteron abundance is lowered. 
In turn the temperature, Tbbn, at which BBN proper 
commences is also lower because the reactions such as 
D + D — > T + p will not begin to process the deuterium 
until a slightly lower temperature. The Universe is there- 
fore older at the beginning of BBN thus permitting even 
more neutron decay than would be expected from just 
the change in neutron lifetime. For 4 He this effect has 
the same sign as the previous two, as Sm np increases the 
amount of 4 He decreases. In contrast, a lower Tbbn ac- 
tually leads to an increase in D because there will be less 
destruction of this nucleus but the increase is insufficient 
to compensate for the previous two effects. 

The effects of a change in Sm np are best summarized 
in figure (|3J) where we plot the deuterium abundance, 



therefore expect more neutron decay in the period from 
neutron/proton freeze-out until the inception of BBN and 
consequently less 4 He and less D. 



Secondly, the neutron-proton mass difference Sm np 
sets the neutron to proton ratio prior to the cessation 
of the interconversion reactions. From equations l|27(l 
and (|28|1 we can see that increasing Sm np will lead to 
increases in both rates per particle but this is more than 
overwhelmed by the decrease in the equilibrium neutron 
abundance F oc exp(—5m np /T). The increase in Sm np 
therefore also leads to lower primordial D and 4 He abun- 
dances because there are less neutrons around that can 
form these nuclei. We can add some quantitativeness to 
these remarks by following the simplified treatment of 
neutron /proton freeze-out of Bernstein, Brown & Fein- 
berg |2fJ. By ignoring n-decay 3 these authors derive an 
expression for the freeze-out abundance of neutrons Y* 
as 

eventually in order to normalize the interconversion reactions. 




Yd, the (free) neutron/proton ratio, F, and the helium-4 
mass fraction, Y , as function of the photon temperature 
at ASm np /Sm np = and ASm np / Sm np — 0.3. At high 
temperatures, above ~ 1 MeV, the shift in the neutron- 
proton mass difference is seen as the difference in the two 
neutron/proton ratio curves while the change in neutron 
lifetime appears as the different gradients for these curves 
close to T ~ 10 5 . The different amounts of neutrons that 
survive from T np to Tbbn results in the large shift in 
the helium-4 mass fraction. This is primarily due to the 
change in the neutron lifetime but is enhanced by a small 
shift in the onset of BBN proper as witnessed in the small 
shift in the position of the peak deuterium abundance 
seen in the top panel. 



D. Constraining 8m np And Bo 

We can use the discussion in the last two sections to 
show how a degeneracy between Sm np and -Bd in the 
prediction of the hclium-4 mass fraction occurs and, fur- 
thermore, illustrate why using deuterium is/could be the 
ideal foil. Raising Sm np will lower the neutron/proton 
freeze-out ratio and so, in turn, reduce the helium-4 
mass fraction but if we also increase Bd then this can 
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FIG. 3: The deuterium abundance, Y D (top panel), the (free) 
neutron/proton ratio, F, (bottom panel, dashed lines) and 
the helium-4 mass fraction, Y, (bottom panel, solid lines) as 
functions of the photon temperature when A8m np /Srn„ p = 
and A8m np /Sm np = 0.3 with a baryon/photon ratio is 6.14 x 



be entirely compensated by an earlier inception of BBN 
proper, the point where neutrons become stabilized in- 
side nuclei. At the same time the predicted abundance 
of deuterium is expected to fall as we increase Sm np but 
an increase in Bd also leads to less deuterium . Thus, 
while the prediction for the abundance of deuterium also 
possesses a degeneracy between the two parameters, it 
is orthogonal to that of helium-4. Therefore, using both 
allows us in principle, to gain strong constraints. This 
orthogonality is best illustrated by figure (@J where the 
complementarity of the two nuclei is evident. 

The figure also shows where the approximations be- 
gin to fail: significant deviations begin to occur when 
ABd/Bd ~ 0.7, which, using the binding energy for 
the deuteron of 2.22 MeV, is remarkably early. The ex- 
planation lies in the Q-value of the reactions we have 
used. During the phase of NSE we require the bind- 
ing energies of tritium and hclium-3 in order to compute 
the reverse reactions and for these we used the present 
measured values. For reference the binding energies are 
-Bt = 8.48 MeV and Bs He = 7.72 MeV. The abundances 
of tritium and helium-3 while they are in NSE are so 
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FIG. 4: The iso-abundance contours D/H of deuterium (hor- 
izontal)and iso-mass fraction contours of helium-4 (vertical) 
from our calculation as a function of the fractional change in 
the neutron-proton mass difference Sm np and the deuteron 
binding energy Bd at 77 = 6.14 x 10~ 10 . From top to bottom 
the deuterium contours are: 4 x 10~ 5 , 3 x 10" 5 and 2 x 10~ 5 
while from left to right the helium-4 contours are 0.25, 0.24 
and 0.23. 



small that even if this introduces a considerable error 
the final results will not reflect this fault. The Q-value 
for the reaction D + D <-» n + 3 He is Q — i?3 Hc — 2Bd 
while that for D + D p + T is Q = B T - 2B D . So 
as we increase the deuteron binding energy the Q-values 
for both reactions decrease and at a 70% increase the Q- 
value for D + D •(-> n + 3 He reaches zero. With such a 
low Q-value the inverse reactions are significant and our 
use of tritium and hclium-3 as neutron sinks is no longer 
valid. In our numerical calculations, when the Q-values 
of these reactions became small, we saw a very different 
flow of the nuclei through the reaction network compared 
to standard BBN. If we persist with the increases in Bd 
and make the Q-values negative we enter very dangerous 
territory since our simple rescaling of the cross-sections 
cannot still apply. Endothcrmic reactions are very dif- 
ferent from exothermic at low energy/tempertaure. This 
is not to say that the deuteron binding energy cannot be 
greater than 1.7 Bd during BBN, it is simply a statement 
that we cannot reliably predict the abundances in this do- 
main. But a 70% increase in the deuteron binding energy 
cannot be regarded as a safe upper limit to the permit- 
ted variation in Bd . At 50% the change in the deuteron 
binding energy is ~ 1 MeV. Since we do not know if the 
change in Bt and -B3 Hc is correlated or anti-correlated 
with Bd, if £?t and -Ba He vary with similar magnitude to 
£>d ! the Q-value of D + D <-> n + 3 He may already have 
been forced to zero at a 1 MeV increase in Bd and our 
approximations cannot be used. Therefore we estimate a 
50% increase in Bd as a limit to the permitted variation 
of this parameter, in the absence of further information 
about how the cross sections and binding energies vary 
with pion mass. 

In contrast, it appears we can vary 5m np in the region 
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we want to explore without running into similar difficul- 
ties. There is a lower limit to the variation we can allow: 
when 5m np drops beneath the electron mass, a decrease 
in 5m np of 60%, neutrons cannot decay so that our rescal- 
ing of the reaction rates with the neutron lifetime to cor- 
rect for the finite mass and thermal radiative corrections 
cannot be applied. Again, a neutron-proton mass differ- 
ence smaller than the electron mass during BBN cannot 
be ruled out, it is just that we cannot calculate what 
occurs when this happens. 

We are now in the position where we can begin to 
constrain our two parameters Sm np and Bd by com- 
paring the predicted abundances with observation. The 
primordial abundance D/H of deuterium is taken to be 
D/H = (2.6 ± 0.4) x 10" 5 |H3 while we use the Olive, 
Steigman and Walker [57| value of Y = 0.238 ± 0.005 
for the helium mass fraction Y. The exact primordial 
abundances remain a topic of debate with two, largely 
incompatible, determinations for the helium mass frac- 
tion j^, |63, E3> 03 |65| and excessive scatter in the 
measurements of deuterium [6fl l69l | but these two nu- 
clei still represent the best probes of BBN because the 
other nuclei that could be used, such as 3 He and 7 Li, 
suffer from large uncertainties in the derivation of their 
primordial values. From comparing the observed abun- 
dances of D and 4 He and their associated errors with 
the iso-abundance contours in figure JIJ it is apparent 
that helium-4 will be the chief source of constraints on 
5m np while deuterium will play the same role for Bd- 
The errors for these observations also reflect the level to 
which we must beat down the systematic errors in order 
to avoid contaminating our results with large offsets. For 
helium-4 we have succeeded handsomely since 0.005 rep- 
resents a 2% error on 0.238 while our systematic was at 
1%. For deuterium we have not done so well since the 
observation has an error of 15% and our systematic was 
at 50%. 

So in addition to the observational errors <td and ay we 
must also include into the analysis the systematic errors, 
Sd and Sy, we introduced when we made our approxima- 
tions. The systematic errors differs from the statistical 
observational errors in that they are correlated. The co- 
variance matrix, V, is therefore of the form 

V - \ S D S Y a Y +S Y J- (32 > 

If we knew the exact predictions for Yjy and Y at any 
combination of 5m npj Bd and r\ then we could easily cal- 
culate Sd and Sy by comparing the exact value with that 
achieved with our approximations. Unfortunately this is 
not the case and we can only make this comparison at 
ASm np /Sm np — ABd I Bd = 0. From comparing our re- 
sults with that from a standard BBN code we can repre- 
sent the systematic errors as the product of the fractional 
errors and the predicted values. We take the fractional 
errors to be the same at arbitrary 5m np and Bd- With 
this understanding we can approximate the covariance 
matrix at all 5m np and Bd and calculate a likelihood via 



£= * exp^T 1 *) 03) 

where the vector 5 is 5 — {Yd — Yd,Y — Y} and the 
hat indicates the observed values. If we also allow rj to 
vary then we have three adjustable parameters and only 
two constraints so we choose to fix this quantity at the 
value of r\ = 6.14 x 10" 10 given by the WMAP observa- 
tions [6(| . The CMB, which is essentially 'atomic physics 
meets cosmology', is not expected to show any depen- 
dence upon Aqcd but other fundamental constants rel- 
evant to atomic physics have been constrained by using 
it WA l6Sj . At this fixed value we derive the con- 
straints on 5m np and Bd shown in figure JSJ . The figure 
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FIG. 5: The 95% and 99% \ 2 contours derived from using a 
primordial deuterium abundance of D/H = (2.6 ± 0.4) x 10 -5 
and a primordial helium-4 mass fraction of Y — 0.238 ±0.005. 
The cross, located at ASm„ p /Sm np — 0.07 and ABu/Bu = 
0.29, indicates the position of the best fit point. 

shows that Sm np and Bd are not as well constrained as we 
might hope because of the large systematic uncertainty 
in the prediction for deuterium but nevertheless the fig- 
ure clearly indicates that the primordial abundances are 
only compatible along a narrow band in the 5m np - Bd 
plane. 

It is worthwhile exercise to show how much stronger 
the constraints would be if deuterium were not contam- 
inated in this way. To this end we show in figure © 
the 95% and 99% % 2 contours when we rescale the deu- 
terium down by 50% and helium up by 1% and remove 
the systematic errors from the covariance matrix. In this 
way we can simulate the situation of a complete under- 
standing of how the cross sections and the binding ener- 
gies are related and there were no reason to truncate the 
reaction network. The size of the contours is now deter- 
mined by the errors in the observational values we used 
and we can see that the current deuterium abundance, 
which primarily constrains Bd would permit only a vari- 
ation of 20% in this parameter at 95% confidence while 
the neutron-proton mass difference would be constrained 
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to within 4%, again at 95% confidence. BBN can pro- 
vide meaningful constraints on the extent to which Sm np 
and Bd could differ from their current values at the ear- 
liest epochs of the Universe if we could determine, with 
more reliability, how the nuclear data we need for the 
calculation depends upon either these two parameters or 
on the underlying fundamental constants. Finally, we 
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FIG. 7: The 95% and 99%y 2 contours as in figure © plus 
the limits upon Bd from |2gl |2S|| after changing the variable 
from the pion mass to the neutron-proton mass difference. 
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FIG. 6: The 95% and 99% % 2 contours derived from using a 
primordial deuterium abundance of D/H = (2.6 ±0.4) x 10 
and a primordial helium-4 mass fraction of Y = 0.238 ± 0.005 
after rescaling the deuterium and helium-4 results from our 
numerical calculation by the systematic errors we identified 
from comparison with a standard BBN code. This methodol- 
ogy is expected to produce results that are representative of 
the situation where we do not have to terminate the nuclear 
reaction network at tritium and helium-3. 

show how we can begin to remove portions of the Sm np - 
Bd parameter space by using their relationships with the 
fundamental constants such as Aqcd- We would like to 
use the relationship for the deuteron binding energy as a 
function of the pion mass derived in Beane and Savage 
|28ll29j| and compare it with our constraints. Beane and 
Savage quote their results in terms of the pion mass, but 
they are varying the ratio of the quark mass to Aqcd- 
They show a range of Bjj vs. which is effectively 
a function Bd = Aqcd f(\Zrn q / Aqcd), where the func- 
tion contains unknown, but constrained coefficients. We 
use this function to vary the QCD scale and relate it to 
the neutron-proton mass difference given in equation 1241 
i.e. 

ABd = AAqcd + Af (34) 
Bd Aqcd f 

In figure (JJJ we superimpose their limits on Bd upon 
our results from figure JSJ. Beane and Savage only con- 
sidered pion masses smaller than its current value so 
while the two sets of curves seem to overlap in a small 
area close to A5m np /5m np = ABd/Bd = a much 
larger overlap might be expected if heavier pions were 
considered. The figure shows that even with the large 
systematic uncertainty in deuterium BBN can rule out 
much of the region of increases in 5m np and Bd up to 



ABd/Bd = 0.5 because it is incompatible with the Bd 
- 5m np relationship from Beane and Savage despite the 
large uncertainty also found in that calculation. New 
limits will be found when heavier pions are considered, 
and also when we have sufficient information to explore 
ABd/Bd > 0.5. 

V. CONCLUSIONS 

BBN presents a golden opportunity to study possible 
changes in the fundamental constants of nature, partic- 
ularly those related to the structure of nuclei and their 
interactions, at one of the earliest epochs of the Uni- 
verse. We have examined the impact of variations in 
one of these constants, Aqcd, upon the predictions for 
the primordial deuterium abundance and helium-4 mass 
fraction. A change in Aqcd will manifest itself through 
shifts in the properties of the nuclei such as the neutron- 
proton mass difference and the deuteron binding energy. 
We have shown how these two parameters are crucial for 
determining the predictions for D and 4 He and how we 
can simplify BBN to the extent that it is a function of 
only these two quantities. While the simplification gives 
very good predictions for the mass fraction of 4 He when 
compared to a standard BBN code the results for D were 
offset by 50%. This large systematic error in deuterium 
swamped the statistical error associated with the obser- 
vation allowing the degeneracy between 5m np and Bd in 
the prediction of helium-4 to show through. By simulat- 
ing the case when the systematic error in the prediction 
for deuterium can be removed we found that BBN can 
limit their variation to the 10% level. In order to make 
BBN a better probe of the time variation of Aqcd, we 
need to know in particular, the dependence of the 3 He 
and T (and 4 He) binding energies on the pion mass and 
the dependence of all the binding energies in the cross 
sections. Even without this input, much stronger con- 
straints are obtained when the result of the BBN cal- 
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culation is combined with the results of the Beane and 
Savage calculations for the deuteron binding energy. If 
Kqcd is the only constant that varies with time, and we 
only consider increases in the pion mass, then in order to 
be compatible with the results from Beane and Savage 
the variations of Bd and Sm np are limited to ~ 1%. 
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